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Several relevant thermodynamic observables obtained within the (2+1) flavor 
and spin zero NJL and PNJL models with inclusion of the 't Hooft determinant 
and 8q interactions are compared with lattice-QCD (1QCD) results. In the case 
that a small ratio R = ^F- ~ 3 at the critical end point (CEP) associated with the 
hadron gas to quark-gluon plasma transition is considered, combined with fits to 
the 1QCD data of the trace anomaly [ 1 1, subtracted light quark condensate HI and 
continuum extrapolated data of the light quark chiral condensate [2], a reasonable 
description for the PNJL model is obtained with a strength g\ ~ 5. ..6 x 10 3 
GeV -8 of the 8q interactions. The dependence on the further model parameters 
is discussed. 

PACS numbers: PACS: ll.10.Wx; 11.30.Rd; 11.30.Qc 

In recent years the role of effective chiral Lagrangians has grown as an impor- 
tant indicator of the order and universality class of phase transitions, as well as of 
the nature and location of the related CEP that may occur for the ground state of 
QCD in presence of external parameters, such as finite temperature T, baryonic 
chemical potential fj,B, magnetic field B [3]. In parallel, 1QCD advances at zero 
and moderate chemical potential with masses approaching the physical values of 
the light quarks [2] and pion mass S, strongly indicate at a crossover transition 
from the hadronic to the quark-gluon phase at finite T and hb = 0. Combining 
1QCD and chemical freeze-out data from relativistic heavy-ion collision facilities, 
the location of the CEP is presently conjectured to eventually occur at R = ~ 2 

and J ~ 1, 0,(61. 
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We consider the SU(3) flavor and spin-0 Nambu-Jona-Lasinio model (NJL) 
(71 with inclusion of the U (\)a breaking 't Hooft flavor determinant iTSTl - lITOTl and 
eight quark (8q) interactions [ 1 1 ],[ 12] (of which there exist two types, one of them 
violationg the OZI rule, with strength g{), and extend it to include the Polyakov 
loop (PNJL) lfT3l - ll22l . The 8q have been firstly introduced to stabilize the ef- 
fective potential of the model ifTTI . Their role turned out to be of significant im- 
portance in the behavior of model observables in presence of external parameters 
[23]-[26],[ 18]-[20|. Of particular interest is that the 8q coupling strengths g\ can 
be varied in tune with the Aq interaction strength G without changing the vac- 
uum condensates and low energy meson spectra, except for the er-meson mass m a 
which decreases with increasing g\. Fits to the low lying pseudoscalar and scalar 
meson spectra yield m a ~ 560 MeV for g\ = 6000 GeV 8 and m a ~ 690 MeV 
for gi = 1500 GeV" 8 02). In the fj,,T plane (where (J, = f-) the g\ , G inter- 
play gives rise to a line of CEP, starting from the regime of large ratios R ~ 20 
(NJL) and R ~ 10 (PNJL) in the case of weak 8q coupling g\, to small ratios 
for strong g\. In the first case the chiral condensate is related with spontaneous 
symmetry breaking (SSB) driven by Aq interactions, in the second scenario SSB 
is induced by the Qq 't Hooft strength 021, EH ASH- This continuous set of CEP 
is particular to the 8q extension of the model.However a correlation between m a 
and the location of the CEP is also observed in the (2+1) - flavor quark-meson La- 
grangian, where besides the 't Hooft term, a quartic mesonic contribution is present 
(27), (28), thus bearing a resemblance to the semi-bosonized version of the 8q NJL 
Lagrangian lfl~2~1l. In order to restrict the g\ values one may: i) calculate decays and 
scattering in the vacuum, which are expected to narrow the choice and ii) compare 
with available 1QCD data at finite T and moderate fi. In the present study we try to 
explore the second option. For the PNJL case an extra uncertainty arises due to the 
parameters related with the choice of Polyakov potential U-p. In particular the To 
parameter of [15], [ 16] has a sizeable effect on the transition temperature. First we 
show in Fig. 1 the CEP lines in a (//, T) vs. g\ diagram. The PNJL model (solid 
lines) enhances the effect of pushing R to small values as functios of g± in compar- 
ison with the NJL case (dashed lines).The crossing of the CEP(T) and CEP(fi) 
lines, ( yielding R = 3), is reached for the PNJL at gi ~ 6.4 x 10 3 GeV -8 , for the 
choice To = 190 MeV, whereas it occurs for the NJL only at a much larger value, 
gi ~ 8.4 x 10 3 GeV~ 8 (we remind that with increasing g\ the crossover becomes 
sharper and eventually gives rise to a first order transition at fi = 0, which happens 
at g 1 ~ 9 x 10 3 GeV~ 8 in the NJL case). Changing T , the CEP{T) is shifted 
up (down) with increasing (decreasing) To (see caption of Fig. 1), while CEP(fi) 
remains sensibly unaltered. In Fig. 2 the chiral condensates and dressed Polyakov 
loop for u, s quarks are shown for the NJL as function of T for fi = and with 
varying strength gi (see caption). For gilO -3 = 1; 5; 6.5; 8 GeV~ 8 the transition 
temperatures T t defined at the corresponding inflection points of the curves are 
T t = 192; 163; 147; 135 MeV for the u-condensate, T t = 197; 163; 150; 135 MeV 
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for the u-quark dressed Polyakov loop, T t = 197; 160; 147; 135 MeV at the first 
inflection point of s-quark condensate, T t = 270; 240; 235; 225 MeV at its 2nd 
inlection point, T t = — ; 166; 150; 135 MeV at 1st inflection point of the dressed s- 
quark Polyakov loop and T t = 270; 240; 235; 225 MeV at its 2nd inflection point. 
The 1st set of inflection points in the case of the s-quark condensate and dressed 
Polyakov loop occur due to the gap equations that correlate the u and s variables, 
yielding similar T t for the u and s observables. The 2nd inflection point occurs at 
temperatures T t larger by ~ 80 MeV A similar pattern is observed for the PNJL 
model in Fig. 3, the second inflection points occur at roughly 55 MeV higher T t 
values. Visually these 2nd inflection points can barely be detected, the transition 
is very slow and smooth. This behavior can be traced back to the fact that for large 
T the s-quark constituent quark mass approaches asymptotically its current quark 
mass value, which is much larger than for the u-quarkQ It is a disputable matter 
which temperature should be taken to characterize the transition for the s-quark 
in these observables. A calculation of the chiral and quark number susceptibilities 
associated with the s-quark in the NJL model display only one peak characterizing 
the transition temperature. In Fig. 4 one sees however that in the PNJL case two 
peaks can occur again for the s-quark chiral susceptibility. Fig. 5 (a) shows the 
trace anomaly calculated for g\ = 6000 GeV 8 , for various values of the parame- 
ter To in comparison with lattice data. In Fig. 5 (b) the subtracted condensate A/ s 
is shown for several values for g\, calculated with To = .19 GeV, and compared to 
1QCD. In Fig. 6 the light chiral condensate is compared with 1QCD data extrapo- 
lated to the continuum limit [2] for different values of g\ and with U-p of |[T5l for 
the cases T = .15 GeV (left) and T = .19 GeV (right). 

From these comparisons we conclude (i) that the smaller the ratio R = ^ 
related with the CEP location, the larger the 8q interaction strength g\ must be 
chosen; a sizeable dependence on the To parameter of the Polyakov potentials can 
induce shifts of the order of several tens of MeV in T c (Fig. 1). For R = 3 we get 
9l of the order 6000 GeV" 8 and T c = 158-188 MeV for the range T = 190-270 
MeV. (ii) Besides the 8q strength, the Polyakov loop plays also a substantial role 
in decreasing the ratio R. (iii) The observables calculated at fi = related with the 
light quarks, chiral condensates, traced Polyakov loop and dressed Polyakov loop 
(Fig. 3), chiral and quark number susceptibilities (Fig. 4 (a), (d)), as well as the 
s-quark number susceptibility (Fig. 4 (e)) and Polyakov loop susceptibility (Fig. 
4 (c)) yield a crossover temperature T ~ 179 MeV for g\ = 6000 GeV -8 and 
To = .19 GeV. (iv) Some of the s-quark observables show two possible transition 
temperatures, Fig.2,3,4(b), the first close to the u-quark transition, the second about 
50 MeV higher for the PNJL model, (v) The best fit to the trace anomaly is for 
gi = 6000 GeV" 8 at T = .21GeV (Fig. 5 (a)) and for the observable A Zs we 



We calculate the thermodynamic potential with the prescription of [26|,|22] where we show 
that it leads to the correct large T asymptotic behavior for the quark masses (condensates), 
traced Polyakov loop and number of degrees of freedom. 
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Figure 1 : Pairs (T, fi) of CEP as function of the 8q interaction gi . Positive slope lines (red online) 
show T-dependence, negative slope lines (blue online) show fi dependence. All model parameters 
fixed as in 1261 . except for gi, G. PNJL potential from 1 15]. Intersection (R = 3) of PNJL curves 
(solid lines) occur at Qi = T = 158; 167; 188 MeV) with gi = 6436; 6251; 6127 GeV" 8 for 
T = 190; 210; 270 respectively (shown only for T = 190 MeV). Intersection of NJL curves 
(dashed lines) at fi = T = 117 with gi = 8372 GeV" 8 . 




the dressed Polyakov loop E l as functions of T for NJL; Figure 3: The same as in Fig.l for the PNJL model. In up- 
solid lines for light quarks, dashed for the strange quark. Up per panel <j> (curves growing with T) stand for the Polyakov 
to down curves in upper panel: g x x 10~ 3 = 1; 5; 6.5; 8 loop and gi strength reverts order compared to the chiral 
GeV -8 , corresponding to black, blue, violet, red (color on- condensate, 
line). 



obtain a reasonable fit with gi = 5.. .6 X 10 3 GeV" 8 and T = .19 GeV (Fig. 
5 (b)). (vi) The peak positions and heights of the continuum extrapolated light 
quark chiral susceptibility vary considerably (Fig. 6). This big spread allows to 
accomodate a large range of g\ values, whose peak positions in turn depend also 
on the choice of the To parameter. The value g\ ~ 5 x 10 3 GeV~ 8 is eventually 
the best choice if one takes the height of the peak also into consideration. 
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Figure 4: Susceptibilities for PNJL, Uv of 02) and To = 190MeV. a) light quark chiral suscep- 
tibilites, (b) s quark chiral susceptibilities, (c) Polyakov loop susceptibility, (d) quark number (for 
u quark) and (e) quark number (for s-quark) susceptibilities. The peaks get more pronounced with 
increasing g\. Color code as in Fig. 2. 




Figure 5: Left: trace anomaly with 8q strength gi = 6000 GeV~ 8 , with Uv from 1 15 ] and respec- 
tive parameter To = .19, .21, .23 GeV (which yield peak positions from left to right). The 1QCD 
data is taken from 1 1 1. Right: the 1QCD data for A; s , the subtracted chiral condensate value normal- 
ized to its zero T value, as defined in 1 1 1, compared to PNJL calculations for several gi strengths, 
color code as in Fig. 2. Solid lines: with Uv from [15], dashed lines: with Uv from 1161 , both at 
T = .19 GeV. 




Figure 6: The 1QCD data for the light quark chiral susceptibility X l chi in the continuum limit 
taken from |2|, in comparison with the PNJL model with Uv 1151 at To = .15GeV (left panel) 
and To — .19GeV (right panel) for different g\ strengths (solid lines, narrower peaks correspond to 
increasing g\). 
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